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Abstract 

Sources of long wavelength radiation are naturally described by an effective field theory (EFT) 
which takes the form of a multipole expansion. Its action is given by a derivative expansion 
where higher order terms are suppressed by powers of the ratio of the size of the source over the 
wavelength. In order to determine the Wilson coefficients of the EFT, i.e. the multipole moments, 
one needs the mapping between a linear source term action and the multipole expansion form of 
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the action of the EFT. In this paper we perform the multipole expansion to all orders by Taylor 
expanding the field in the source term and then decomposing the action into symmetric trace free 
tensors which form irreducible representations of the rotation group. We work at the level of the 
action, and we obtain the action to all orders in the multipole expansion and the exact expressions 
for the multipole moments for a scalar field, electromagnetism and linearized gravity. Our results 
for the latter two cases are manifestly gauge invariant. We also give expressions for the energy flux 
and the (gauge dependent) radiation field to all orders in the multipole expansion. The results for 
linearized gravity are a component of the EFT framework NRGR and will greatly simplify future 
calculations of gravitational wave observables in the radiation sector of NRGR. 
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I. INTRODUCTION 



The multipole expansion is a standard tool in physics which students usually encounter 
first when studying electrostatics [l| , where it is employed to derive approximate expressions 
for the electric field at a distance R from a source of size a provided that a <^ R . In this case 
it is an expansion in a/R and is employed at the level of the solution of a static equation 
of motion, where its form can be derived from Taylor expanding the Green's function of 
the Laplacian. In time dependent situations, the multipole expansion does not correspond 
to an expansion in a/R any more. For the radiation field far away from the source in the 
regime R S> a, A, all multipoles yield a leading contribution to the far zone radiation field 
falling off as 1/ R. For the leading 1/R piece of the radiation field, the multipole expansion 
is organized as an expansion in a/ A which may be truncated for long wavelengths A ^> a. 

One important feature of the multipole expansion is that the multipole moments are 
organized in irreducible representations of the rotation group SO (3). This has several im- 
portant advantages. It makes calculations simpler and more transparent, and it ensures 
the absence of mixing of multipole moments for example in the energy flux in a linear the- 
ory (so that we can truly speak of dipole radiation to all orders, otherwise there would be 
for example dipole-trace-of-octupole radiation at higher orders). Performing the multipole 
expansion then includes a decomposition into multipole moments which are in irreducible 
representations of SO (3), where we use symmetric trace free (STF) tensors. 

While the traditional field of application of the multipole expansion has been electro- 
magnetism, it is used in general relativity and in particular in gravitational wave physics 
extensively (J 4]. In both cases however, it is commonly applied at the level of the solution 
of the equations of motion when computing the fields. In this paper, we want to study the 
multipole expansion from a different angle, where we consider the multipole expansion at the 
level of the action. This is required in many nonrelativistic effective field theories (EFTs) 
in order to obtain a uniform power counting {5, 6]. 

EFTs have become an indispensable tool of modern theoretical physics to systematically 
study systems with multiple separate scales, see |7|] for a review. Their construction is 
based on the physical scales in the problem and the underlying symmetries, which together 
provide the form of the action. It is important to have a uniform power counting for all of 
the ingredients of the theory, i.e. at the level of the action, so that one can systematically 
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calculate to a given order without having to guess which terms or Feynman diagrams need 
to be included. The coefficients in front of the terms (or operators in quantum language) 
in the action are called Wilson coefficients and need to be determined from matching. In 
weakly coupled theories where the underlying "full theory" is known the matching can be 
oerformed analytically. Examples for such EFTs include NRQED, NRQCD [8j and NRGR 
6] which were constructed to describe nonrelativistic bound states. 

In these theories, the size of the bound state, i.e. the size of the radiation source, is 
much smaller than the wavelength of emitted radiation. In order to disentangle the different 
scales and to achieve a uniform power counting, it is necessary to introduce separate modes 
of the fields which describe physical effects on different scales and which have different 
kinematic properties. Potential modes for example yield the leading binding dynamics of 
the source while radiation (or ultrasoft) modes describe physical on-shell radiation. In order 
to have a uniform power counting of the action written in terms of the different modes, 
it becomes necessary to Taylor expand the radiation modes describing the physics at the 
longest distances in the action around a point, which lies within the source . This is 
the origin of the multipole expansion in these nonrelativistic EFTs, where it is required at 
the level of the action and plays an important role in the matching to an effective radiation 
theory. 
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The EFT framework NRGR [f| (see for a pedagogical review and 
for nonrelativistic gravitational bound states such as compact binary systems has become an 
important systematic tool to perform post-Newtonian (PN) computations for gravitational 
wave templates which are crucial for direct detection experiments such as LIGO/VIRGO 
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12| or eLISA 



13] . It is an alternative to the more traditional methods to perform PN 
calculations, see [14J for a review. 

The EFT construction of NRGR to describe compact binary systems is based on the 
hierarchy of length scales during the inspiral phase, where the size of the compact objects is 
much smaller than their orbital separation which in turn is much smaller than the wavelength 
of the emitted gravitational waves. The first effective description used to describe the physics 
at distances between the finite size scale and the orbital scale is constructed by replacing 
the compact objects with point particle worldlines, where their finite size effects can be 
included with higher order terms in the point particle action. The second stage of the EFT 
is constructed for scales between the orbital and the wavelength scales, where the potential 
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modes which are responsible for physical effects at the orbital scale are removed from the 
theory. The resulting theory has a conservative part and a radiation sector in terms of the 
radiation modes describing gravitational waves. 

The radiation sector of NRGR has been further developed and explored beyond the lead- 



ing order in [15j, where the general form of the action for the radiation sector was constructed. 
Its form is determined by the underlying symmetries, reparameterization invariance and dif- 
feomorphism invariance, and is applicable to arbitrary gravitational wave sources in the long 
wavelength approximation. This action of the effective long wavelength radiation theory is 
in the form of a multipole expansion and is a derivative expansion where higher order terms 
are suppressed by powers of the ratio of the size of the source over the wavelength. The Wil- 
son coefficients of the action, the multipole moments, are not determined by the symmetries 
and need to be fixed through a matching calculation. 

This matching onto the effective radiation theory is performed perturbatively using Feyn- 
man diagrams. In particular, we calculate in the "full theory", i.e. in the effective theory 
which still has potential modes as degrees of freedom, Feynman diagrams with one external 
radiation mode and arbitrarily many potential modes 1 being exchanged between the two 
point particles. The result of these Feynman diagrams can be written as an effective linear 
source action S = —l/(2mpi) J d^xT^h^. Here, T^ v is the stress-energy pseudo-tensor 
which includes for example the gravitational binding energy of the system, and we define it 
in terms of the Feynman diagrams. This action is certainly not in multipole expansion form 
with a derivative expansion. In fact, to obtain a uniform power counting, it is required that 
the radiation field in this action is Taylor expanded. This Taylor expansion around a single 
point forces the resulting radiation theory to be formulated in terms of just one worldline 
coupled to gravity. So in order to match onto the effective radiation theory, we need to 
first Taylor expand this action, write it in terms of manifestly gauge invariant operators and 
finally express the coefficients in irreducible representations of SO (3). This then gives the 
multipole moments in terms of moments of T^, which in turn are computed using Feynman 
diagrams. 

In [jjj] this matching was first done to NLO (for the quadrupole) and resulted in repro- 
ducing the energy flux to 1PN. In doing this, the decomposition of the coefficients of the 



1 Which diagrams need to be included to a given order is fixed by the power counting of the EFT 
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Taylor expansion of the radiation fields into irreducible multipole moments needed to be 



performed, which at that order was a straightforward exercise. At higher orders [16| . this 
decomposition becomes increasingly cumbersome, and the purpose of this paper is to solve 
it to all orders. Our results will streamline the matching needed in higher order calculations 
of gravitational wave observables 2 in the effective field theory framework NRGR. 

In this paper we perform the multipole expansion at the level of the action first for the 
simple case of a scalar field, then we consider electromagnetism and finally we study general 
relativity with a linearized source term. While our results for gravity are our main result, 
we demonstrate our methods used in multipole expanding the action first in the simpler 
cases. For all three theories studied, we multipole expand the action to all orders and give 
the exact multipole moments. We also compute the energy flux or radiation power and the 
radiation field far away from the source in terms of the multipole moments. 

The multipole expansion at the level of the solution of the equations of motion has been 
performed for the same three cases we consider here by Damour and Iyer in Q . Their results 
are given for general sources, but taking the long wavelength limit of their expressions, we 
find complete agreement with our results for the multipole moments. 

Notation: We work with a mostly negative metric signature (+, — , — , — ). Greek let- 
ters denote Lorentz indices ranging from ... 3 and lower-case latin indices denote spatial 
indices from 1 ... 3. Upper-case indices denote a set of spatial indices where we use the multi- 
index notation introduced by Blanchet and Damour [18]. For example, x N = x kl . . .x kn or 
x ijN-2 _ x i x i x ki _ _ _ x kn ~ 2 . Parenthesis around indices denote symmetrization, for exam- 
ple A^B j C k) = | (A i B j C k + A i B k C j + A j B { C k + A j B k C i + A k B i C j + A k B j C { ) . We use 
natural units c = H = 1. 



2 NRGR and NRGR-inspircd worldline EFTs have produced many interesting results in gravitational wave 
physics and beyond, see 
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II. SCALAR FIELD 



A. Multipole expansion of the action 



We first study a scalar field 4> coupled linearly to a source J. Its action reads 



S = J d 4 x y-d^(j> + J<j)j (1) 
and its equation of motion is 

U(j)=J. (2) 

Now we consider the configuration where the spatial variation of the field outside the source 
is much larger than the size of the source J, i.e. if spatial derivatives scale as di<p ~ ^<f) 
and if the size of the source is a, we are considering the situation a < A. For radiation 
this means we are working in a long wavelength approximation. In this case we can Taylor 
expand the field <fi in the source term in the action Source around a point within the source, 
and we choose our coordinates such that the point we expand around is the origin x = 0. 
We plug the Taylor expansion 

CO -. CO ^ 

W> X ) = E n\ ^ • • • **" {8kl ■ ■ ■ 8kJ) (t ' 0) = E "j xN (*> °) ( 3 ) 

n=0 ' ra=0 

into the source term in the action 

//* n OO -4 n OO -4 

d i xJ<P = dt d 3 xJ(t,x)^-x iV (^)(t,0)= / dt^2~M N d N <P (4) 
J J n=0 U ' J n=0 U ' 

where we defined the moments M N (t) = j d 3 xJ(t,x)x N . Note that all M N are already 
symmetric in their indices k% . . . k n . In order to bring the source action into the form of a 
multipole expansion we need to decompose the moments M N in irreducible representations 
of the rotation group SO (3) for which we use symmetric trace free (STF) tensors. 

The tensors M N can be expressed in terms of STF tensors starting from the formula for 
arbitrary symmetric tensors S N 0, 
[»] 

qN _ qN | n\(2n — 2p — 1)!! x ( fclfca s:k 2 p-ih 2 p qk 2p +i...k n )a 1 a 1 ...apa p /e\ 

b - b STF + 2^ („_ 2 p)!(2n-l)!!(2p)!! " ( ) 

where [|] denotes the largest integer < n/2. On the RHS of Eq. (J5J) the tensors of lower 
rank, gtop+i-fcnaiai-apa^ s ^iH are no t trace free in their free indices k 2p +\ . . . k n . Therefore, 



we will use Eq. (jSJ) recursively to make all tensors trace free, where we mean that all free 
indices in {k\, . . . , k n } which are not on a 5 kikj should become trace free. This results in 

[f] 

p=0 

where the STF prescription only applies to the uncontracted indices /c2 P +i • ■ • k n and where 
the coefficients are 

Jn) = n\(2n - 4p + 1)!! 

p (2p)\\(n - 2p)\(2n - 2p + 1)W 1 ' 

With this the source action becomes 



*EEn- d3 * J ° r2Px sTF P (v 2 ) p ^_ 2P 
= / * g H § (2^t|W / (V2) ' ^ (8) 

where we define r = |x|. Now we may use the equation of motion outside the source □</> = 
to convert the contracted spatial derivatives to time derivatives, and in turn integrate by 
parts to let them act on the moments rather than on the field. This yields the source action 
in multipole expanded form 
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u source 



r 00 1 

J dt^-l^ (9) 



oc 



^ (2p)\\ (2£ + 2p+ 1)!! J * stf v j 



1=0 

with the multipole moments given by 

(2£ + iy. 

(2rAV. 

p=0 

The normalization of the multipole moments II is chosen such that for p = their expression 
reads 

X L = J rf 3 xJxf TF , (11) 

where we note that the static case of a time independent source is given by only the p = 
components of the multipoles. 

B. Energy flux for scalar radiation 

We compute the power or energy flux emitted in scalar waves to all orders in the multipole 
expansion from Eq. (J5J) using the method employed in {9], Q . For the amputated amplitude 



of (j) emission from the £th multipole moment we have 

iA® = itf^h* (12) 
where = (|k|,k) is the outgoing momentum of the emitted on-shell scalar 0. We then 



compute the energy flux as in [9|, 



E = ± [ ^-\A\ 2 = — ! dQ * dk k2 y y Xk l+i l L (k)T l \k) n^ni (13) 
2T J (2tt) 3 1 1 2T J (2tt) 3 ^ ^ £\£\ 



£=0 i =0 

where we defined rik = k/|k|. The angular integral is performed with 

if p odd 

<mn iJ ={ (14) 
^pjy §( klk2 . . . § k p-i- k p) if p even 

which results in the condition £ = £ since the multipole moments are trace free and for 
£ I some 5 will contract two indices of the same multipole moment. Moreover, when using 
Eq. (fl4j) . we need to omit the combinations of indices on the #'s which vanish because the 
multipole moments are trace free, or in other words, we only need to keep the products of 

T f* 

o's where each 5 contracts one index on I with one index on X . That is accounted for by 
a factor of (£\) 2 2 e / (2£)\ = £\/{2£ - 1)!! inserted into Eq. ^ and contracting the multipole 
moments with each other. We find 

1 />oo °° 1 

E - / dk V 1 f, 2ii+1) T L T L * 

E -^Tj dk 2^£ K 2£ + l)ll k 11 

= £ 477i(27TT)!! ( (dt^ lL ) ) (15) 

where the first line is the expression for the energy flux in momentum space and the second 
line is the energy flux in coordinate space. In the latter, <> denotes time averaging. 



C. Radiation field far from the source 



At a distance R ^> A ^> a only the 1/R component of the field is relevant and we may 
neglect the components which fall off faster than 1/R. The field <f)(t, x) can be computed 
from the one-point function, i.e. the 0-emission amplitude with a retarded propagator for 
the (f> leg. The use of the retarded propagator in diagrammatic field theory calculations can 
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be derived in the in-in formalism [20J. For classical calculations of the field (i.e. the one- 
point function), one can show that the use of retarded propagators is the only modification 
the in-in formalism implies, even if nonlinearities and multiple couplings to a source are 
considered. Let us show how to isolate the leading 1/R piece of the multipole expanded 
radiation field. For that purpose it is convenient to write the multipole expanded action of 
Eq. © as 

A£^xV)^)^---^/\y), ( 16 ) 

1=0 ' y y 

and we require the retarded propagator which in d=4 spacetime dimensions reads in coor- 
dinate space 

-iO(t f -U) 5(f _ ^ {) ^ (17) 



47r|xj — x vi 

Then the field or the one-point function is given by 

««■ x > = < l d> yt\ I'WW^ • • • ^T^f a(t " «° " Ix - y|) . (18) 

Note that the spatial derivatives can either act on the delta function 5(t — y° — |x — y|) or 
on the factor l/|x — y|. For one spatial derivative acting on these two possibilities we have 

dy k \x-y\ 

—S( X ° - y° - | X - y|) = d ( X ° - l X ~ y ^ - 5( x ° - yO _ | x - y|) 

dy k " ' dy k d(x° — |x — y|) 

Mx -y -\x-y\)) . (20) 



(x- 


y) k 




x — 


y| 3 




d(x° 


— X 


-y|) 




Qyk 




(x- 


y) k I 


' d 


x - 


-y| v 


v dy 



We see that if we integrated these expressions over y using the S 3 (y) in the expression of the 
field in Eq. ( IT8l) . only the second one where the spatial derivative acts on the delta function 
would give a term proportional to l/|x|. When further derivatives act, they can also act on 
the vectors (x — y) fc in Eqs. (|T9| [20]) . These however will give traces which vanish since the 
multipole moments are trace free. After integrating over y, all derivatives acting on powers 
of |x — y| will yield terms which fall off faster than l/|x|. Thus, all spatial derivatives V y 
have to act on the delta function, where they are effectively converted 3 to — n ^o where 



3 This corresponds to computing the amputated one <f> emission amplitude on-shell in momentum space, 
since this corresponds to setting k = fc°n. 
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n = x/|x| . If we then employ this in the expression for the field we find for the radiation 
field to all orders in the multipole expansion 

1 00 1 

<P(t, x) = 4^hd ^J\ nL ^ JL )^) ( 21 ) 

where we integrated by parts to let all time derivatives act on the multipole moments and 
where the moments and their time derivatives are evaluated at retarded time t ret = t— |x|. We 
can also use the radiation field of Eq. (l2Tj) to compute the energy flux using E — f dVt\x\ 2 cf 2 
which confirms our result of Eq. (fT5l) . 



III. ELECTROMAGNETISM 

A. Multipole Expansion of the action 

The action for the electromagnetic field coupled to a source J^ 1 reads 

S = J d A x {-^F^ - J»A^ (22) 

where F^ u = d ii A v — d v A^ and where the current in the source term is conserved, 

<V M = . (23) 

The electric and magnetic fields in terms of the potentials are E = — VA° — A and B = V x A 
respectively and obey the usual Maxwell equations which are given outside the source by 

V-E = 0, V-B = 0, VxE = -B, VxB = E. (24) 
Plugging the Taylor expansion of the field into the source term of the action S source , it reads 

Source = - I dt f rf 3 xJ^(t, x) J2 -7 a^A, • (25) 

J J n=0 U ' 

As the first step, we will express Source in terms of gauge invariant operators. This means 
that, aside from the monopole term arising from the /1 = and n = piece, all other 
components of Eq. (125]) will be given in terms of the electric and magnetic fields and their 
spatial derivatives. For this purpose, we need to consider the A and the A 1 components 
separately. The former are simply written as 

SiL» = -J dtQ d 3 xJ^A° + I d 3 xJ°x N ^d N ^{-V kn A°). (26) 
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and we see that, aside from the first term coupling the total charge of the source to the 
electrostatic potential, the A source action includes couplings to the A piece of the electric 
field E. 

Now the A % piece of the source action, 



J*±h(J- 

J n =0 w 



d 3 ^J i x N ) d N A 



(27) 



needs to be rewritten such that it provides couplings to the magnetic field as well as the 
A 1 components of the electric field which were missing in the A part in Eq. (126 p . This is 
achieved by decomposing the source action for A 1 into two parts, each with definite symmetry 
properties in the spatial indices {i, k\ . . . k n }. The first part is totally symmetric in all indices 
and the second one is antisymmetrized in one pair of indices (i,ki) and then symmetrized 
in {ki . . . k n } . These symmetrizations are conveniently done using Young symmetrizers 2l| 
denoted as Young tableaux. In terms of these, the couplings in Eq. (127)) are expressed as 



cf x J l x 



N 



(n + 1)! 



% 


h 




k n 



n 



(n + 1)! 



% 


k 2 




k n 


h 





+ k-perms^ 



(28) 



where 



fci|.. Jfcnl = n\ J rf 3 x [fx N + + k-perms)] 



(29) 



i 


k-2 




k n 







+ k-perms = n\ / d 3 ^J i x N - (n - 1)! / rf 3 x (j^ar** 2 "^ + k-perms) , (30) 



and where "+ k-perms" means that all other combinations of the indices {ki . . . k n } which 
differ from the original are to be added. For example, it means that in Eq. ()28[) there are n 
Young diagrams to be added, each with a different index kj in the box in the second row. 
With this decomposition the action becomes 



n=0 



1 



n\ \ (n + 1)! 



k\. . .k n + 



(n + 1)! 



i 


k 2 




k n 


h 





•)) 



■ k-ponus j j ()nA\ 
(31) 



For the totally symmetric piece we now use the conservation of J M from which follows that 

1 



J d 3 xj°x iN = -J d 3 xd a J a x iN = J rf 3 x [fx N + (J k ^x ik2 - kn + k-perms)] = -^ | j {hj . \k 



(32) 
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i 


k 2 






h 





and the part with the antisymmetrization can be written as 

+ k-permsj = (n - 1)! J rf 3 x (e*^*/!*-^ + k-perms) . (33) 
Using these expressions in Eqs. fl32|) and fl33l) . the A* part of the source action gives 

=j dt f:^l d 3 xJ°x N d N ^ (-A^) -J dtjr J d 3 x (J x x)*-^"^^*- 



(34) 



where we integrated by parts to let the time derivative act on the field A 1 . Now we can 
combine this the source action for A adding Eqs. ( 126]) and ( 134]) and obtain 

Source = - J dt (^J d 3 xJ ^) A° + Jdtjr^Jj J d 3 xJ°x N d N ^E kn 

„ oo „ 

- J dt E (^T)j J ^ x ( J x x ) fc " ^dN-iB^ , (35) 

which is in the desired form. Clearly, the source action is now in a manifestly gauge invariant 
form, however the coefficients of the operators are not yet in irreducible representations of 
SO (3), i.e. they are not yet in STF form. Bringing the action in a form with coefficients 
in STF for is performed by first writing all coefficients as symmetric tensors and secondly 
taking out the traces using the formula of Eq. (J6]). 

The first term in Eq. ( 135]) is trivially in STF form, whereas the coefficients of the other 
terms in the first line of Eq. ( 135]) are already symmetric in all indices {k\ . . . k n }. Thus, we 
only need to bring the terms in the second line of Eq. fl35|) , 



„ oo „ 



xx)*"/" 1 ^, (36) 



in the form of symmetric tensors. For this purpose we rewrite the decomposition in terms 
of Young tableaux of Eq. ff28l) in the convenient form 

K i X N ~ 1 = [K l X N ~ l ] S + - ^ic e abc K a x bk 2 ...k n ^ + k _ perms ) (37) 

where we made use of Eq. f )33|) . Using this in the action gives 

s 



i n + 1)! 

n=l v ' 



J x x) k " x N ~ l 



d N - X B h 
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n=2 



n 



n — 1 



n + 1)! n 



(fx ((J x x) x x)"- 1 2 <9 jV _ 2 (-V x Bp- 1 



^ (n + l)! 

n=l 



:j x x) fc - /- 1 



A. ; ?l. 



- / dtS^ 

/ ^ (n + 2 

J n=l K 
n. oo 



d X 9f J • X X 



N 



n=l 



n + 2 ! 



(i 3 x d,J kn r 2 x N - 1 



d N -xE k 



(38) 



where we have used the Maxwell equations outside the source to replace the curl of B, and 
we integrated by parts to let the resulting time derivative act on the moments rather than 
on the fields. The last term in Eq. (I3"8"|) is still not symmetric, and in order to bring Source 
into a form expressed only in terms of symmetric tensors, the symmetrization procedure 
of Eq. fl371) has to be applied infinitely many times where we use the Maxwell equations 
in vacuum to replace any curls of E or B and integrate all time derivatives by part to 
act on the moments. This is straightforward since the resulting structures are reoccurring 

(2) 

cyclically, and can we derive the expression for Source in terms of symmetric tensors. With 
this decomposition, the total source action in terms of symmetric tensors is 

s 



SfL 



(It 



d 3 xJ° A 



J n=l 



d 3 *J°x N 



d N -iE k 



n=l q=0 

„ oo oo 

+ /*££ 

J n =l q=0 

oo oo 

<*££ 
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n + 2q + 2)\ 

n 

'n + 2q + 2)\ 



d 3 x d? q+1 J ■ xx w r 29 



d 3 x d} q+1 J k -x N - l r w 



d N -!E k 



n 



[n + 2q + l)\ 



d N -iB h 



(39) 



d^8 2 t q (J x x) fc " x N ~ l r 2q 

III •■• I M / 

n=l q=0 

Now the final step is to use Eq. ([6]) to write the symmetric tensors in terms of STF 
tensors. The first term in Eq. (1391 coupling the static potential to the total charge, 



S Q 



dt I / d 3 xJ° A 



(40) 



is trivially symmetric and trace free. For the other parts, we can essentially repeat the 
procedure of the scalar case. However, we need to omit the traces which are zero on-shell, 
i.e. which give terms in the action proportional to V • E = or V • B = 0. Moreover, we 
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need to account for different structures of the moments which can arise from traces of the 
moments. Let us start with the non-trivial part of line 1 of Eq. ( 139]) . 

s 



J n=l 



d^rx 



(41) 



When using Eq. ([6]) for the trace free decomposition of we need to count how many of 
the symmetrizations of indices in Eq. fl6]) give terms with V • E. For a tensor with n — 2p 
free indices on Sstf and 2p indices on <5's on the RHS of Eq. (jHJ) there are (") possible 
combinations of indices which are on the 5's (not counting their possible orders). Out of 
these possibilities, there are (""-J possibilities of indices on the 5's which include the index 
k n and which is the index on E kn in the action. Thus we need to subtract these possibilities 
when we set V • E = 0. This gives the piece of the action 



*EE-„. 

n=l p=0 

oo _ oo 

*EfE 



Cp \2p) \2p-l) 



dWr 2 ' X ^ p (V 2 )" d N -2 F -lE k ~ 



(2£+ 1)!! 



+ 2p (2p)\\ (2£ + 2p + 



— J d 3 x d?J°r 2 *>xi TF d L ^E k * , (42) 



where we used the wave equation DE = to convert V 2 E to E and then integrated by parts 
to let the time derivatives act on the moments rather than on the field. 

The trace free decomposition and the counting for the piece of the source action in the 
second line of Eq. fl39|) works exactly like the one for S El mce because in both cases all free 
indices of the moments are on x . We obtain 



nE 2 



oo oo 



*EE 

n=l q=0 
oo oo / p 

*EE(E 

1=1 p=0 \j=0 



n + 2q + 2)\ 



ci 3 x df q+1 J ■ x/r 2 « 



d N ^E h 



(£ + 2p + 2)\ 



rf 3 x (df p+1 3 ■ x) xl TF r 2p d L - l E k * . 

(43) 



However when we take traces in the STF decompositions of the remaining terms of Eq. 
( 139]) . we need to count how many times a factor of r 2 , J • x, or (J x x) • x = occurs. The 
components of the source action in the third line of Eq. fl39l) can be expressed as 



cE 3 

*— ' ar\i it 



oo oo 



*EE 

n=l 5=0 



n 



n + 2q + 2)\ 



dN-\E k 
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oo oo L 2 J 
n=l q=0 p=0 



n 



n) n-2p 



n + 2g + 2)! p n 



n 



jk n _ 2px N-2P~l r 2{q+p)+2 + j _ xx Af- 

n 



-2P r 2(q+p) \ 



STF 



x oo / p (^+2j) ' 



£=1 p=0 \j=0 



+ 2j i (£ + 2p + 2)! 



d 3 x a 2p+ W'V 2 ^ 



STF 



„ oo oo / p r> • (£+2j) 

+ / *££(£ 



= i p =o \i=o 



+ 2j / (£ + 2p + 2)! 



rf 3 x (<9 2p+1 J • x) x^pr^d^E^ . 



(44) 



The counting of terms from the STF decomposition for the fourth line of Eq. (|39|) works 
essentially the same way as for the third line, only that now (J x x) • x = and the only 
possible structure to arise from contracted indices is r 2 . We find 



s?± 



n=l q=0 



11 



[n + 2q + 1)! 



kn x N ~ l r 2q 



oo oo / p (£+2j) \ 



«EE E 



1=1 p=Q 



-£ + 2j J (£ + 2p+l) 



J d 3 x<9 2(? (J X X 

J rf 3 x d 2 t p ( J X X 



STF 



dL-lB k < 



(45) 



Finally, we can add Eqs. gOJ), (|42]), (T43J), ()44]) and (T45J) to find our final result for the 
multipole expanded electromagnetic source action 



Sr. 



| dt (| d 3 xJ°) + | dffjil^E* + J dtf^j^J L d L ^ 



(46) 



with the exact expressions for the electric and magnetic multipole moments 



'' = £ 

p=0 
oo 

+ £ 



(2£+l)!! £ 



(2p)\\(2e + 2p + £ + 2p 
(2£+l)\\ 



d 3 xd 2p J°r 2p x L 



n STF 



p=0 



(2p)!!(2£ + 2p + 1)!! (£ + !)(£ + 2p + 2) 



y d 3 xr 2 ^ 2p+1 (JV-V 2 - J ■ xx L ) 



jL = Y1 



(2£+l)H 



p=0 



(2p)\\{2£ + 2p + 1)! 



/ rf 3 x 9 2p (x x j; 



STF 



-i STF 

(47) 
(48) 
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We note that we used Eq. (1A2j) to eliminate the sums over j present in Eqs. (1431 - [45]) . 
and we chose a normalization 4 such that the electric moments in the static limit where 
p = are 1^ = 

x stf an d the magnetic moments for p — are given by Jq 

J d 3 TL (x X J) 



X 



L-l 



STF 



. Due to the conservation of the electromagnetic current d^J^ = 0, 
our expression for the electric multipole moments in Eq. f)47p is not unique, since we can 
add or subtract terms which vanish due to current conservation. In particular, current 
conservation implies 

s 



d 3 xd+J°r 2p x L 



J d 3 xr 2p - 2 (iJ^x^r 2 + 2pJ ■ xx L ) 



(49) 



which can be used to write p > pieces of I L in the first line of Eq. f T4"T|) in terms of the 
form of the components seen in the second line and vice versa. 



B. Radiated power 

We calculate the energy flux or radiated power in the same way as above for the scalar 
field in Sec. Ill Bl In order to compute the total energy flux however, we need to sum the 
amputated photon emission amplitude squared over all polarizations h, and one needs to 
choose a gauge to have an explicit expression for the polarization sum. Since it is slightly 
simpler we will employ transverse gauge or Coulomb gauge, where the polarization vectors 
satisfy e°(k, h) = and k ■ e(k, h) = and where the polarization sum is 

e*(k, h)e j *(k, h) = S ij - W /k 2 . (50) 

h 

In this gauge there is no A emission, and the amputated emission amplitude is 

iA » = 1 E ~ir lL{ W) ^We^k, h)+iJ2 77XnT ji (l k D ^" W^V*( k , h) . (51) 
e=i l - i=i ^ + >■ 

Squaring the amplitude, summing over helicities with Eq. (|50|) and performing the angular 
integration with the help of Eq. fTl4l) yields the energy flux 



1 f°° 
E = —— / dk 
^ 2 TJ 



OO 



V — - k 2(e+1) \l L (k)\ 2 + V - k 2{l+l) I J L (k)\ 

^ee\(2£ + i)v. 1 1 j| + A-(£ + i)!(2£ + i)!! /c |J l/cj| 



4 Note that our normalization differs from some of the usual conventions used in electrodynamics (see e.g. 
[1]) because we adapt the conventional normalization for multipole moments in gravitational wave physics. 
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v e+1 !( de+1 A 2 \iT e l( de+1 A 2 \ (52) 

j^4n££\{2£+l)\\ \V^ +1 / / ~ 4vr(£ + 1)!(2£ + 1)!! \ V^ +1 ) / K ' 

It is straightforward to check that if we use the covariant Feynman gauge instead where we 
need to keep A , we obtain the same result for the energy flux. This is of course expected 
because our multipole moments are gauge independent and the energy flux is an observable 
and cannot depend on the gauge. 



C. Electromagnetic radiation field 

In calculating the radiation field we will work in Feynman gauge where the retarded 
photon propagator is given by 



4/t|xj — Xj| 



■8(tf - ti - |x/ - Xj|) . (53) 



Using the one-point function we can then calculate the radiation fields A°(t, x) and A l (t, x) 
as discussed in more detail for the scalar field in Sec. Ill CI We obtain for the multipole 
expanded radiation field in Feynman gauge 

oo 



1=1 



A'(f,x) 



47T X 



OO . CO n 



(55) 



IV. LINEARIZED GRAVITY AND NRGR 



For gravity our conventions are m Pl = 1/(32ttG) for the normalization of the Planck mass 
relative to Newton's constant, the Riemann tensor is R^vafi = d a Y^ — 9pT^ a + T\aFvp ~ 
r^r^ a , the Ricci tensor is R^ u = R a IWi v l and we remind the reader that we work with 
a mostly negative metric signature. When we expand the metric around flat Minkowski 
spacetime, we normalize the gravitational field as g^ v = r]^ u + h^/mpi. 

The action we want to study is given by the standard Einstein- Hilbert term plus a linear 
source term. It reads 

S = -2m 2 Pl [ d 4 x^R - — *— I d A xT» v h^ , (56) 
J 2m P i J 
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where d^T^ = 0, and the vacuum equations of motion are 

iV = 0. (57) 

The linear gravitational source T^ v may be either the standard stress-energy tensor of 
matter or it may be the stress-energy pseudo-tensor which includes some gravitational ef- 
fects such as the gravitational binding energy of a bound state. The latter is the case in 



the effective field theory framework NRGR 



6|, where different kinematic regions or modes 



of the gravitational field are treated separately by writing h^ u = h^ u + H^ u . The radi- 
ation modes h^ v describe on-shell gravitational waves whereas the potential modes H^ v are 
responsible for the conservative binding dynamics of a bound state such as a binary system. 

The characteristic length scale for the potential modes is the size of the source a, whereas 
for the radiation modes it is the wavelength A ^> a. Therefore if we consider physical effects 
at distances much larger than the size of the source a, the potential modes are not needed as 
explicit degrees of freedom in the theory. Integrating out the potential modes, i.e. solving 
for them and replacing them by their solutions in the action, we obtain an effective radiation 
theory only in terms of the radiation modes. The effective action of this effective radiation 
theory is expanded in powers of the radiation field h^, 

S cS = S + S X + 0{h 2 ) (58) 

and is obtained perturbatively using Feynman diagrams. The action So of zeroth order in 
hfxv, obtained from Feynman diagrams with no radiation modes and with any number of 
potential modes exchanged between the two wordlines describing the two binary constitu- 
tents, gives the action describing the conservative dynamics of the system. The first order 
couplings in Si are obtained from Feynman diagrams with any number of potential modes 
being exchanged and one external radiation mode and can be written as a linear source 
action 5 . Thus when we use the term linearized gravity in this paper, what we mean is that 
the source term in the action is only coupling linearly to the radiation field. In other words, 
we treat the coupling of the effective source T^ v to the long wavelength modes linearly, 
while the interactions of short- distance modes (which are implicit in T^ v ) can be included 
to arbitrary order. 



5 Terms in the effective action with coupling of quadratic and higher order in the radiation field become 
relevant for higher order calculations, in post-Newtonian applications they start at order v 5 or 2.5PN. 
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A. Multipole expansion of the action 



As in the previous sections, we Taylor expand the gravitational field h^ v in the source 
action Source around a point within the source choosing our coordinates such that the point 
we expand around is the origin x = 0. The source term in the action then is 

5 sour ce ="-!— / d A xT^h^ = --^— [ dt V / d 3 xT^(t, X )x" (d N h, v ) (t, 0) . 

ATripi J ^rripi J ^ n\ J 

(59) 



As in the electromagnetic case, our first step will be to express this Taylor expanded source 
action in terms of manifestly gauge invariant operators. In general relativity, the gauge 
invariant quantities are the Riemann tensor, the Ricci tensor and the Ricci scalar and their 
covariant derivatives. Due the vacuum equations of motion = R = however, only 
structures uwo.ving the Rieuranu tensor yield physica. tenus which can contribute* 

observable. Terms with Ricci tensor or scalar can be removed via field redefinitions [6|, 12 2| ; 
in EFT jargon such terms are called redundant operators. Now when writing the Riemann 
tensor in terms of temporal and spatial indices, there are only three distinct combinations 
of indices: Roioj, Roijk and Rijki- In linearized gravity however, the purely spatial Riemann 
tensor is on-shell (i.e. upon use of the vacuum equations of motion) 

Rijkl — —t'ijmCklpRomOp j (60) 

so that we may replace it by -Ro?oj- The couplings to the two remaining components of the 
Riemann tensor, Roioj and Roijk, are conveniently organized in terms of the electric and 
magnetic components of the Riemann tensor 6 . They are given by 

= RoMj = (d djh 0i + d dih 0j — didjh 00 — d^h^) (61) 

E{j ~z£imnRojmn 7} ^imn (d()3 n hj m -\- <9j<9 m /l()n) > (62) 
Z ZTflpi 

where the last expressions are expanded to linear order in h^ v . Note that we gave the 
expressions for and in a system with four- velocity = (1,0,0,0), which makes 
sense for our application to the multipole expansion where we Taylor expand around the 



6 On-shell, where R^, the Riemann tensor R^ vp<J and the Weyl tensor C^ pa used in previous work |15| 
coincide. 
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point x = which implies the source is described by a worldline x^ = (t, 0) whose four- 
velocity is v" = (1,0,0,0). 

The electric and magnetic components of the Riemann tensor in vacuum have the prop- 
erties Eij = Eji, En = 0, Bij = Bji and Bu = 0, i.e. they are symmetric and trace free. 
Furthermore, derivatives acting on the electric and magnetic components of the Riemann 
tensor in vacuum give the equations of motion which are in complete correspondence to 
Maxwell's equations, 

d{Eij 0, diBij 0, £i mn d m Ej n Bij, £imndmBj n Eij , (63) 
which imply the wave equations 

DEij = UBij = . (64) 

Having reviewed these properties of the electric and magnetic components of the Riemann 
tensor, we set out to write the Taylor expanded action in terms of E^ and B^. Similar to 
the electromagnetic case where the coupling to the total (conserved) charge could not be 
expressed as a coupling to E or B, we also expect that the couplings to the leading multipoles 
in gravity cannot to be expressed in terms of E^ and B^. 

In order to express the Taylor expanded action mostly in terms of E^ and B^, we have 
to consider the couplings to hoo, ha and hij separately. First we write the hoo part of the 
source action as 

S^Z CC = / dt [(/ d^T^ h 00 + (/ rf 3 xT°v) d k h 00 

+ / *£ ^ (/ d 3 xT 00 (t,x)x N ^J d N _ 2 (-L- (-c^.^oo)) • (65) 

where the first two terms are the couplings to the total energy and to the center of mass 
coordinate (times the total energy) respectively. The h 0i part reads 



qhoi / jf 

^source ' <" 

mpi 

oo 



d 3 xT 0i )hoi+[ / d 3 xT 0l x k d k h 0i 



(66) 



n=2 

where the first term is the coupling to the total linear momentum and the second term 
contains the total angular momentum. All terms besides the first one need to be decomposed 
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to forms of definite symmetry using Young tableaux, where we may use our results from the 
E&M calculation for the vector potential in Eqs. f)28l- [30]) with J 1 replaced by T°\ For the 
totally symmetric piece we then use the conservation of T^ v from which follows that 



d 3 x.f 00 x iN = - I d 3 xd a T Ua x 



: J d 3 x [T 0i x N + (T 0kl x ik2 - kn + k-perms)] = ^ | j \kx\. . \k 



and for the part with the antisymmetrization we use 



i 






k n 


h 





™ + k-perms 



in 



1)! / c/ 3 x (e ifelc e afec T°V fc2 - fc " + k-perms) 



(67) 



(68) 



Note that Eqs. (!67|) and ( 168]) and also are identical to Eqs. ( 132|) and ( 133|) in the decompo- 
sition for the electromagnetic case, with (J°, J a ) replaced by (T 00 ,T 0a ). 
With these the part of the action becomes 

s%L* = J dt (J d ^ T ° l ) hoi+ {\J rf3x ~ T ° v ) } \ ( dkh ° l ~ dlhok) 



+ j dt fl ^ (/ dW ) 9 "-z {dodk^ho^ + dod^ho^) 



n=2 
oo 



n=2 



2n 



n + l)\ 



d ^ e k n ba T 0a x b kl ...k n ^ | q 



W-2 



2m P i 



(efc„cd5fc n _i<9 c ^od) 



For the part of the source action, 



ntiij 
w source 



the decomposition in Young tableaux reads 

1 



[ dt lt-i f d 3 ^x N d N h ij , 
2mpi J n J 

^ l J n=0 J 



(69) 



(70) 



d 3 xT tJ x 



N 



+ 



(n + 2)! 

n — 1 
2(ra+ 1)! 



k\ . . .k n + 



71+ 1 

fn + 2)! 





.7 


A-2 




fen 







+ k-perms 



i 


j 


h 




k n 




k 2 





+ k-perms 1 . 



(71) 



The first Young tableau is given by 



i 


j 


fci 




k n 



2n! / (< 3 x[ ^ + ((W + r^-) +k -Pe n ns) + (r—,«- + k - P e rms )] 
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= n\ J d 3 xf 00 x ijN } (72) 
where we used the conservation of T^ u to derive the last equality, and the second one is 
+ k-perms 



i 


3 


k 2 




k n 


h 





(n-1)! / d 3 x 



in 



2nT ij x N + (n - 2) ((T^ar** -1 + T^x*" -1 ) + k-perms) 
- 4 (T k - k ^ x ijN ~ 2 + k-perms) 
1)! J d 3 x [ e ik " c e abc {T aj x bN - 1 + (T akn ~ 1 x bjN - 2 + k-perms) } + (i j)] + k-perms . 

(73) 

In the last line of Eq. (l73p . the first +k-perms inside the parenthesis means that each of the 
indices ki which are not on the epsilon tensor, i.e. which are inside the round parenthesis, 
is to be used once as an index on T aki inside the round parenthesis. The second +k-perms 
in Eq. (IT3I) means that each ki is to be used once as an index on the epsilon tensor. The 
third Young tableau is 



i 


3 


k 3 




k n 


h 


k 2 





+ k-perms 



2 (n-2)! / cfx 



n(n - l)T lj x N -{n-1) ((T ik "x jN ~ l + T jkn x N ^) + k-perms) 
+ 2 (T k " k ^ x ijN ~ 2 + k-perms) 



2 (n-2) 



When we plug these into the hij part of the source action we obtain 



-perms 



(74) 



J n=2 n ' •* 



2m P i 



(-d 2 h kn _ lkn ) 



0a„6JV-l a 
X CaT-2 



n=2 
oo 



2mpi 



(efe„cd^0^d^cfc n _i) 



+ / ^Ef^-TV / ^l 7 ™*" ' + T k -^ 2 x N - 2 -2T k " a x N - la }d N . 2 E kn _ lkn 

J n =2 ^ ^ >' *' 

(75) 



where we have used 



J d 3 xf 0a x bN = J d 3 x [T ab x N + {T ak -x bN - 1 + k-perms)] 



(76) 
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to arrive at the expression in the second line of Eq. (!75j) . and we used the epsilon tensors 
from Eq. f!74|) to write dk n _ 1 dk n hij as the purely spatial Riemann tensor, which in turn is 
written as the electric components of the Riemann tensor via Eq. (I6T?|) . 
Adding all components of the source action, the total becomes 



C _ ocons _i_ orad 

^source ^source ' ^source 



(77) 



with 



neons 
^source 



2m P i 
1 



dt 



2m P i 
1 



dt 



2m 



pi 

OO 



xT uu J h 00 + 2^J cfxT Ui J hoi 
dt [Mh 00 + 2P l h 0i + MX%h 00 + Ve^kfyhok] 



T°V) ) - (d t h ok - d k h 0i ) 



(78) 



C1 



rad 



+ 



■' n=2 

oo 

dtj2 

n=2 

oo 



n — 1 
(n + 1)! 

2n 



d 3 x [r/ + T fcn_1 



k n aN—la~] o p 

J CiV-2-^fc n _ifc, l 



n=2 



n + 1)! 



d 3 xe * n6 a r 0a a;6 JV-l ^ 



(79) 



We note that S s c ™ s ce , 



which is not written in terms of E^- and B^, couples to several quantities 
which are conserved due to the conservation of T^ u : The electric parity mass monopole 
M = J d 3 xT 00 is the total energy (ADM mass) of the source 7 , P J = J <i 3 xT°* is the total 
linear momentum (which interestingly does not have a consistent classification as a multipole 
moment since some of its aspects are characteristic of a monopole such as its coupling 
without derivative, but it is a vector like a dipole), and the magnetic parity current dipole 
L 1 ^ = — J rf 3 x (T *x J — T°ix l ) is the total angular momentum (which is related to the usual 
angular momentum 3- vector via L % i = e l i k ~L k ). The only coupling in Source which is not 



00 



X . 



conserved is the mass dipole, related to the center of mass position X* = 1/M J d 3 xT 
However, d a T^ u = implies MX* = P l and X* = 0. Thus, none of these quantities can 
radiate since the energy flux would be proportional to the square of the first time derivatives 
of M and P* or the second time derivatives of the dipoles. Moreover, it is straightforward 
to show that under coordinate transformations S°°™ e only changes by total time derivatives 



7 Note the change in notation from [1 51 ] where the mass monopole was denoted m. 
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which leave the dynamics unchanged. It usually is convenient to work in the center of mass 
frame where X* = P* = 0. 

The couplings in the second part of the source action S™£ Tce are all able to emit radiation, 
and we now set out to decompose Source i n terms of STF tensors. As done in the previous 
two sections, we first write all structures as symmetric tensors. For this purpose we write 



Source — $0 + S^' + S^' + S 2 (80) 

S » = J dt E ^ / rf3x { T °° + ^T™}x N d N _ 2 E kn _ lkn (81) 
'° = " / dt £ / d 3 xrS^ d N ^E kn _ lK (82) 

J n=2 ^ ' )' J 

Sf' ° = - J dt g J d 3 xe k ^T°°x bN -i d N . 2 B kn _ lkn (83) 

S * = J dt lt (^Tlj! / ^ x r^SV^ d N „ 2 E kn _ lkn . (84) 

The subscripts on the terms of the action denote the number of pairs of indices which can 
be antisymmetrized without necessarily vanishing in the action. We see that So is already 
symmetrized. 

The strategy to write the action in terms of symmetric tensors is the following: First, 
we perform the (partial) symmetrization of S 2 reducing it to terms which have at most 1 
pair of indices which may give a non-vanishing contribution to the action when antisym- 
metrized. Then we perform the symmetrization of all terms where one pair of indices can 
be antisymmetrized similarly to what we did for the electromagnetic case. 

For the reduction of of S 2 we use the decomposition via Young tableaux from Eqs. ( 1711 - 
( 1741) which, aside from the last line in Eq. ( 172|) where we used the moment relations following 
from dfjT^ = 0, hold in the same form when additional factors of r 2 are inserted in the 
integrals. Otherwise, the procedure works analogous to the electromagnetic case where we 
decomposed Eq. ( )36|) by using the decomposition with Young tableaux infinitely many 
times. In doing so, we also use the moment relation 

J d 3 x [T ab rV + (T ak "r 2p x bN - 1 + k-perms)] = J d 3 xf 0a r 2p x bN - 2p J d^T ac r 2 ^x hcN 

(85) 



with 
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to rewrite the second Young tableau's expression in a more convenient form. After some 
quite lengthy manipulations, we obtain 

s 



So 



n — 1 



n + 2q + 1)1 
n — 1 



n=2 q=0 

oo oo 

+ I dtYY , 

I L^L^ < n + 2 q + l)\ 

n=2 g=l v ' 

oo oo 



d N - 2 E kn _ lkn 



d 2g / d 3 x.T aa r 2q x 



N 



2(71 — 1] 



(n + 2g + l 

n=2 g=l v ' 



dN-2Ek n ^k n 



d 2q / dPvT^i&x*"- 1 \ 8m-9.E, 



N-2^k n -!k n 



oo oo 



2n 



(86) 

Plugging this into S™^ Tce we have 



<7 



rad 
source 



n=2 
oo oo 



+ /<*££ 

^ n =2 <j=0 
„ oo oo 

+ /*££ 



to"/ 

n — 1 
{n + 2q+ 1)! 

n — 1 



dN~2E kn _ 1 k n 



1 s 



A 29 / rf 3 xT fc„_ 1 /c„ r 2 g+ 2 a;J V-2 



<9jV-2-Efc„_ifc„ 



(n + 2g + l)! 
2(n - 1) 



ft 29 / d 3 *T aa r 2q x 



N 



^Af-2-E'fc n _ifc n 



n=2 q=0 

oo oo 

n + 2g + l 

n=2 <j=0 v ' 
oo oo 2 

rft ££ (n + 2o)!(n + 1) 

n=2 g=0 v H > v y 

"rf 3 xe fc " fea {a^r 2 ^- 1 - 2g5 29 " 1 T ac r 29 - 2 x bc7V - 1 } d^B^u, 

(87) 

where in the sums over q in the last three terms the q = components were provided 
from 5*o, Sf' and S'f' respectively. Only the last two expressions in Eq. flBTj) need to be 
further symmetrized, where the procedure works analogously to the E&M calculation. For 
the fourth line of Eq. ( 1H71) the symmetrization procedure yields 

2(n - 1) 



oo oo 



n=2 <j=0 



ra + 2g + 1)! 
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oo oo oo 
n=2 q=0 s=0 



„ oo oo oo 

-/*£££ 

J n =2 q=0 s=0 



2n(n - 1) 



n + 2g + 2s + 3)!(n + 2s + 2) 



s a; a6JV 



(i 3 x Q^ +2s+2 T ab r 2q+2 

2n(n - 1) 
n + 2g + 2s + l)!(n + 2s) 



<9Af-2-Efc n _i/ 



1 s 



oo oo oo 



+ /*£££ 

7 n=2 <j=0 s=0 



and for the last two lines of Eq. flHTj) it gives 

/oo oo „ 
rftVV, ^ r 



d 3 x ^^T^r 2 ^ 2 ^^- 1 d N ^E kn _ lkn 

2n(n - 1) 
n + 2g + 2s + 2)!(n + 2s + 1) 



<i 3 x d 2 i +2s+1 e k ^ ba x ac r 2q+2s x bcN ~ 1 
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(89) 



so that we obtain an expression for S™^ Tce entirely in terms of symmetric tensors. Its form 
can be simplified by removing the two infinite sums over q and s in Sf and Sf using variable 
redefinitions q — > q + s which yield one infinite sum over q and a finite sum over s. The 
latter can all be performed using Eqs. (1A5|) and (IA6I) and we find the final result for Source 
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in terms of symmetric tensors, 
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r„, + $>„ + 2^! / 4 

n=2 q=0 

The final step in the multipole expansion of the action is to take out the traces so that all 
coefficients are in terms of STF tensors. For this purpose, let us define 4 separate components 
of Eq. ( )90l) . The first piece Si comprises the first four lines in Eq. (1901) . the second piece 
S 2 are lines five and six, the third piece S3 is line seven and the fourth piece S4 is given by 
lines eight and nine. Each of these four pieces is written in terms of STF tensors using Eq. 
OH]), but the counting of resulting trace structures differs. 

For the first piece S\ all contracted indices in the STF prescription will give powers of 
r 2 , so that the only counting left to do is to omit terms which vanish by En = diEij = 0. 
This is similar to what we did for the electromagnetic case (where the V • E = components 
had to be omitted, see Eq. f|42|) and the discussion above). Here this is accounted for by 
inserting a factor of 

\2p) \2p-l) \2p-l) 



(91) 



into the STF decomposition of Eq. ([6]) and only keeping the traces which don't yield En or 
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diEij. We obtain 
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For the remaining pieces S2, S3 and S4, we also need to count which of the different possible 
structures of traces in the moments occur how many times. For example, in the T okn part of 
S2, the traces when applying the STF decomposition of Eq. ([6]) can yield either r 2 or T 0a x a , 
and we need to know how many times each of them occurs. It is a straightforward counting 
exercise, and for the parts of the action S2 and £4 which have only one of the free indices 
hi . . . k n not on a vector x, the counting is the same as in Eq. (jUJ). Moreover, we note that 
the factor in Eq. fl9lT) is universal for all terms and is to be used multiplicatively with the 
factors which account for the different traced structures in the moments. We obtain for the 
second part of the action 
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STF 



the third one becomes 
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and the fourth part reads 
S A = - 
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Adding these and simplifying we find for S™^ rce in terms of STF tensors 
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While the above Eq. (|96|) is in the desired form in terms of STF moments, it is still much too 
complicated. In order to simplify and to get rid of the sums over j, we first use the moment 
relations resulting from d tl T f11 ' = to simplify the expressions for the electric moments. In 
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particular, we first use 
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in order to replace the components with T ke ~ lke , and subsequently we apply 



d 3 xT ake r 2p x aL - 



d 3 xf 0a r 2p x aL 

d 3 xT aa r 2p x L 



2p 



d 3 xT ab r 2p ~ 2 x abL 



(98) 



and 
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to replace all structures with T akg and T oke . 

After these replacements, we can perform all remaining sums over j using Eqs. (1A2 
[4]), and arrive at our final result for the multipole expanded source action 
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with the exact expressions for the multipole moments 
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The above form of the multipole moments I L and J L is not unique and can be modified 
using moment relations derived from d^T^ u = 0. See the appendix [B] for the power counting 
rules which tell us which terms are required at a given order in PN calculations using NRGR. 



B. Energy flux in gravitational waves 



In order to calculate the energy flux or radiated power, we follow the same procedure ex- 
plained above in the scalar and E&M sections. We here work in linearized gravity neglecting 



nonlinearities, see 



15] for calculations in the EFT framework of the leading effect of non- 



linearities in the gravitational wave energy flux. We use a physical gauge with polarization 
tensors satisfying e° M (k, h) = e ll (k, h) = kV J (k, h) = and where the helicity sum is [jjj] 
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(107) 

In this gauge there is only emission, and the amputated on-shell graviton emission 
amplitude is 
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(108) 



Upon squaring of the amplitude, summing over helicities using Eq. (11071) and evaluating 
the angular integrals with Eq. (1141) . we find the energy flux 



-tJ dk 2^^-i)m+m 1 m 
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which is of course gauge invariant and reproduces the well-known result of [4]. 
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C. Gravitational waveform 



The gravitational radiation field far away from the source is computed in harmonic gauge 
where the retarded propagator in d = 4 spacetime dimensions is 



■S(t f -U 



*l) 



(110) 



471-lx/ - Xj| 

with P^ap = \ + r]^r] ua — rj^ v rj a p). The physical waveform is transverse traceless 

and is obtained from the harmonic gauge radiation field via contraction with the transverse- 
traceless projection operator. Conforming to the usual dimensionless normalization, we 
define the transverse traceless waveform to be 



mpi 
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with the transverse traceless projection operator given by 



Att„(a;) = (6 lk - n t n k ) (8 jt - n,n ; ) - - (6 tJ - n,n,) (5 kl - n fc n,) . (112) 



We find the waveform to all orders in the multipole expansion in linearized gravity to be 
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(113) 



where we point out that our overall sign differs from the (otherwise equivalent) expressions 
for hj? in [y| because we use a different metric signature. Moreover, we can also derive the 



energy flux of Eq. (11091) from the waveform in Eq. f II 1 3 j) using 
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D. Going beyond a linearized source term 



Our result for the multipole expanded action for linearized gravity in Eq. (UOOp with the 
moments in Eqs. ( 11011 - I106p can simply be extended beyond the linear regime using the 
diffeomorphism invariance of the action. In NRGR this is guaranteed to hold in the action 
of the effective radiation theory by the use of the background field method when integrating 
out the potential modes |6|. 

Operationally, covariantizing the action corresponds to replacing all Riemann tensors 
with their expressions to all orders in the gravitational field h^ u and replacing all partial 
derivatives by covariant derivatives in Eq. fllOUl) . as well as extending the couplings in 
^source beyond the linear re gim e. We then have the form of the action written down based 
on symmetry arguments in [lSj, which in the center of mass frame is 

Source = -M J dT ~\ J dx"Ve ijk col k 

/oo r 00 2£ 

t * r 12ji lLVl -* E *-i*i - / dr E?J-W jL VL-2B ke _ lki (115) 
e=2 ' t=2 ^ '' 

where we can also replace dr = dt^/goo and dx^ = dtd^ when we parameterize the worldline 

in terms of coordinate time t. The Wilson coefficients of this action with couplings beyond 

the linear regime in Eq. f 1 1 1 5 j) are still given by Eqs. (11011 - I106j) . the same multipole 

moments as the ones in the linear multipole expanded action in Eq. (UOOp . 

However, we note that new couplings quadratic in h^ u can be written down which are 

not included in Eq. (II 15ft . for examples terms proportional to EijE k i- Their coefficients are 

not fixed by our work. Moreover, we expect that there exist some relationships between 

the one-graviton-emission amplitudes used to match at the linear level and multi-graviton- 

emission amplitudes, since they must conspire to combine in such a way that they provide 

the non-linear pieces in Eq. All 5j) . 



V. CONCLUSIONS 



We have performed the multipole expansion of a linear source term action in the long 
wavelength approximation for three theories of increasing complexity, namely a scalar field, 
electromagnetism and general relativity. Throughout, we worked at the level of the action 
which we expressed in terms of manifestly gauge invariant operators in electromagnetism 
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and general relativity. We provided the exact results for the multipole moments, and we gave 
the emitted radiation power and the radiation field far away from the source to all orders 
in the multipole expansion. Our results for the multipole moments agree in the appropriate 
long wavelength limit with the results in 2]. 

Our results for linearized 8 gravity are an important component of the effective field theory 
framework NRGR. Our exact expressions for the multipole moments in Eqs. ( 11011 - 1 106[) 
together with the expressions for the power emitted in gravitational waves in Eq. ( I109p and 
the waveform in Eq. (I113P will greatly simplify future calculation of precision gravitational 
wave observables in the effective field theory framework. From now on, all that is left to do 
in the matching to the effective radiation theory is to compute the Feynman diagrams which 
provide the moments of T^ v in the expressions for the multipole moments to a desired order. 
See the appendix |B] for the power counting rules which tell us which terms are required at 
a given order. 

While we have worked with a linear source action in gravity, our results hold more gen- 
erally for the covariantized form of our multipole expanded action in Eq. ( 11 15ft by diffeo- 
morphism invariance. However, at quadratic order in the gravitational field, new invariant 
operators can be written down, and this work does not fix their coefficeints. The multipole 
expansion to all orders in the nonlinearities has been solved at the level of the solution to 
the equations of motion in the more traditional approach to compute gravitational wave 
signatures by Blanchet in j^J], where the introduction of four additional sets of multipole 
moments was required. It would be an interesting future direction to further investigate the 
structure of the source action and the multipole expansion beyond the linear regime and its 
relation to the results of 0. 
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See the discussion in Sec. IIVI for what precisely we mean by linearized. 
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Appendix A: Useful formulas for sums 



Recalling that 
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as denned above in of Eq. (GO), we have the sums 
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Appendix B: Power counting for radiation observables in NRGR 

The multipole moments in Eqs. (11011 - II 06[) together with the energy flux in Eq. (j!09p 
and the waveform in Eq. ( 11131) greatly simplify future higher order calculation in the EFT 
framework NRGR. Here we briefly outline which multipole moments and need to be included 
for a calculation at a given order, and which of the terms in the infinite expressions for the 
multipole moments are required. 

Generally, all time derivatives in the multipole moments or in the expressions for the 
power and the waveform scale asv/a where v is the typical three- velocity within the source. 
All factors of x in the moments scale as the size of the source a (which is the orbital separation 
for binaries). Therefore, the next order in p in the sums in the multipole moments in Eqs. 
(11051) and (11061) . which always is accompanied by an extra factor of d 2 r 2 , is suppressed by a 
factor v 2 . The only ingredient we are missing in order to be able to power count everything 
are the moments. Their scaling is not universal and changes for example if we consider 



systems with spin. See Tabled 16J for their scalings where we defined K^ v = j d?xT^ v x L . 
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e s 
ma v 


£ 6 

ma v 



TABLE I: Scalings for leading terms in moments of T^ u for various orders in spin [16]. These 
scalings are valid for all moments I > 2 and persist if some indices are contracted. 



For practical purposes it is convenient to simply power count factors of v with respect 
to the leading expression for a multipole moment or an observable. Let us discuss a quick 
example for a system where we can neglect spin: If we want for example the quadrupole 
moment to 2PN, i.e. at order v 4 beyond its leading expression, we need the T 00 component 
in Eq. fll05p to second order in p, the T aa and T 0a components in Eq. f ll05p to first order 
in p and the T ab component in Eq. f ll05p to zeroth order in p. Recall that each order 
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